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Based on the photon-exciton Hamiltonian a microscopic theory of the Casimir problem for di- 
electrics is developed. Using well-known many-body techniques we derive a perturbation expansion 
for the energy which is free from divergences. In the continuum limit we turn off the interaction 
at a distance smaller than a cut-off distance a to keep the energy finite. We will show that the 
macroscopic theory of the Casimir effect with hard boundary conditions is not well defined because 
it ignores the finite distance between the atoms, hence is including infinite self-energy contribu- 
tions. Nevertheless for disconnected bodies the latter do not contribute to the force between the 
bodies. The Lorentz-Lorenz relation for the dielectric constant that enters the force is deduced in 
our microscopic theory without further assumptions. 

The photon Green's function can be calculated from a Dyson type integral equation. The geometry 
of the problem only enters in this equation through the region of integration which is equal to the 
region occupied by the dielectric. The integral equation can be solved exactly for various plain and 
spherical geometries without using boundary conditions. This clearly shows that the Casimir force 
for dielectrics is due to the forces between the atoms. 

Convergence of the perturbation expansion and the metallic limit are discussed. We conclude 
that for any dielectric function the transverse electric (TE) mode does not contribute to the zero- 
frequency term of the Casimir force. 



I. INTRODUCTION 

More than 50 years ago Casimir published two completely different papers [1] [2] on what was later called the 
Casimir effect: One [1] was a cumbersome microscopic theory, the other [2] a very elegant macroscopic one. From the 
latter one may get the impression that the problem is rather simple, namely, quantization of the free electromagnetic 
field in a region with boundaries. However, the results of such macroscopic calculations (the so-called mode summation 
method) are notoriously divergent ([3] and the many references therein). Although by suitable regularizations one may 
get reasonable results, the ultraviolet divergences cannot be eliminated in general. This shows that the macroscopic 
problem (with hard boundary conditions) is not well-posed. A second weak point of the macroscopic method is that 
it completely ignores the fact that the real Casimir force is the result of electromagnetic interaction between neutral 
atoms and not a free problem. This was already the point of view of Casimir and Polder in their first paper [1]. 
It is clear from the above considerations that a real understanding of the Casimir problem can only come from a 
microscopic theory. 

In the following we develop such a theory for dielectrics and use it to solve various Casimir problems explicitly. 
The basis is the photon-exciton Hamiltonian given in the next section. The interaction consists of the non-relativistic 
dipole coupling to the radiation field and the dipole-dipole approximation of the Coulomb interaction. In Section III 
we derive the perturbation expansion for the energy E for arbitrarily spaced atoms. This energy is finite since infinite 
self-energy contributions of the dipoles are excluded. 

When taking the limit to a continuous dipole distribution one has to be earful because of the dipole-dipole approxi- 
mation. We turn off the interaction at an average distance a between the dipoles. This was also done by Marachevsky 
[6], leading to a finite energy for a dilute ball. To extract the classical Casimir energy from E we have to introduce 
a further ultraviolet cut-off e~ 5k . E s splits then into two parts: E s = E D — Eq. In the limit 5 — > on\yE s remains 
finite, whereas E D and Eq are divergent. In the limit a — > also E is divergent. E D is for S, a — > formally equivalent 
to the usual Casimir energy Ed as given by Candelas [8]. It is always infinite. 

To calculate the force between separated bodies in regions A and B, one has to use the free energy F. It splits 
into three finite parts: F — Fa + Fb + Fab- Fa and Fb contain the free energies of the two isolated bodies. They 
are finite since the finite distance between the atoms is taken into account and self-energies are excluded. Only the 
interaction energy Fab depends on the distance of the bodies and contributes to the force. In the limit a — > only 
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this energy remains finite. For the polarization we obtain the Lorentz-Lorenz relation without further assumptions. 
The force between disconnected bodies can be calculated by differentiating F with respect to the distance or by a 
contour integral over a region containing one body. Only the stress field outside the bodies contributes to the force 
and is finite also when the cut-off distance goes to 0. 

In Section IV the Green's function of the total electromagnetic field for a continuous dipole distribution field is 
discussed. It is given by a Dyson type integral equation. In this equation the many-body problem is separated from 
the geometry. The latter only enters through the region of integration which agrees with the region of space occupied 
by the dielectric. We identify a class of geometries which allow exact solution of the integral equation. As applications 
we consider in Section V and VI plain geometries and we re-derive the Lifshitz formula for the Casimir force between 
two plates [7] . In Section VII we solve the integral equation for the Green's function of the electromagnetic field for 
a dielectric ball. 

In Section VIII the metallic limit e — > oo is considered. Within our theory because of the Lorentz-Lorenz relation 
this means to put the atomic polarization oeo(u) = 3 for u = 0. It is not clear whether that the perturbation expansion 
converges. We'll show convergence in flat geometries just for u = in the TM mode and for small frequencies u > 
in the TE mode. The TE mode does not contribute in the zero-frequency case u = to the Casimir force independent 
of the model adopted for the dielectric function. This result cannot be deduced from the Lifshitz formula for the 
plasma model where lim„^o e(iu)u 2 ^ 0. The zero- frequency behavior plays an important role in the context of the 
temperature dependence of the Casimir force and is widely discussed in recent papers [11-17]. 

II. THE MODEL 

We consider N atoms in interaction with the quantized radiation field. The dynamical degrees of freedom of the 
atoms are their internal excitations, therefore, our system is described by the following photon - exciton Hamiltonian: 

H = H% + H% h + H u (1) 
H% = Y,E n b+{x)b n (x), (2) 



(: E{xf : + : B(xf : ) 
I d 3 xka+(k)a x (k) (3) 



A=l,2 



* =E[-^W) + ^:^) 2 : 

X 

-\ J2 ^x-x'). (4) 

x^x' 

The atoms sit at places Xj , but for simplicity of the notation we omit the index j and write the sum over all atoms 
as Y]f b^(x) and b n (x) are emission and absorption operators for an atom in state n at the place x. The interaction 
Hi consists of the usual dipole coupling of non-relativistic radiation theory plus the dipole - dipole interaction <!>. We 
assume that the distance between the atoms is big compared with the extension of the atoms (Bohr radius). Then $ 
is given by 

3 

§{x-x') = Qi(x)qj(x')^i,j(x- x') (5) 
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where q(x) denotes the dipole operator of an atom at place x. It can be expressed by one-particle matrix elements of 
the position operator, similarly to the momentum operator p(x): 

P(x) = ^2{fO,p(Pn)bn{x) + {tp n ,P<Po)b+(x). (7) 

n 

The commutation relations for the emission and absorption operators are 

[b m (x),b+(x" )] = S mn S S x', (8) 

and zero otherwise. Sss' is actually a Kronecker delta, it is =1 if a;' = x and otherwise. The quantized electromagnetic 
field is given in the radiation gauge by 



(9) 



where £\(k) are transversal polarization vectors: 

(&(£),&(£)) = <*av, (10) 
(ix(k),k) = o. (11) 

The photon emission and absorption operators satisfy 

[a x (k),a+(k')]=6 X yS(k-k') (12) 

and zero otherwise. 

III. ENERGY AND THE STRESS TENSOR 

We use the Matsubara formalism to calculate the expectation value for the energy for finite temperature. The 
expectation value of Hamiltonian is given by: 

E = (H)= T ^~!ff = -9,mTr( e -^) (13) 

With 

S = e-f 3H (14) 

we obtain 

E = -^m ^-g^ -^lnTr^o) 

= -d fj \n(S) a + (H a ) a . (15) 

The term (H )o is just giving the ground state energy of the free system. It is the sum of the excitation energies of the 
atoms and the black-body radiation energy. The later is for T — > 0, since we use the normal ordering in the definition 
of the free photon Hamiltonian (3). The normal ordering is necessary for the Hamiltonian to be well-defined. The 
(-ffo)o-term will be ignored in the following. The expectation value of the S'-matrix is given by perturbation theory. 

00 f) 

(5)0 = 1- y)(-l) n / dT 1 ...dr n (T T H 1 (T 1 )...H 1 (r n ))o (16) 

n=l J 



where H\{t) is the Hciscnberg representation for the interaction Hamiltonian Hi with imaginary time r = —it. It is 
well known that in ln(5)o only connected graphs contribute [4, 5]. We have 



00 00 ^ 

ln(S)o = ^ ^- ^ a (iw m ,£i).D9 i2 (iw m ,£i -£ 2 ) 

m=-oon=2 xi...x n 

. . .a (iw m ,x n )D° nil (iui m ,x n - xi) (17) 
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where the last sum goes over the positions x% of the dipoles with x\ ^ X2 • • • for not including the divergent self-energy 
contributions, is the atomic polarization given by 

aUj(T,x-x') = -(T T qi(T,x)qj(0,x')) = a (T 7 x)S i: jS xx > (18) 

for isotropic atoms. Since a obeys the KMS periodicity ao(r) — a n (r + 0), it can be expand into a Fourier series 

a {T,x) = - a Q (i^m,x)e' tuJmT , (19) 

m— — 00 

where 

2irm 

u m = — (20) 
are the Matsubara frequencies. The Fourier coefficients in (19) are given by 

f 

&o(w,x) = / drao(T, x)e WT , (21) 
Jo 



using the short-hand notation 

.2-Kvn 



W = VjJ m = I- 



(22) 



The Green's function D^(u),x — y) in (17) is the Fourier transformed free Green's function of the total electric field 
given by 

D9 j (T,x-y)=d^(T T A i (r,x)A j (0,y}) + ^ ij (x-y} (23) 

In momentum space we get 



MM) 



k% kj 

js> 2 ' Jz'i k j lu S j j ~ k j 



,2 



= -^ + ^2 E (24) 

A=l,2 



We get the energy in (15) by differentiating (17) with respect to and using dpuj m — —uj m /0. Hence 

^ 00 00 ^ 

m— — 00 n— 2 x\---x n 

... a (ioj m ,x n )Di nil (ioj m ,x n - £1). (25) 
To calculate the Casimir force with finite temperature we need the free energy given by 



= ~p S E^ 2 ao(«w m ,fi)D? li2 (iw ro ,£i - f 2 ) . . . 

m— — 00 n— 2 xi-.-Xn 

... 6io(iUm,x n )Di nil (iLj m ,x n - xi). (26) 

In the limit T — > (/3 — > 00) the two energies £7 and F are the same. Indeed the summation over the frequencies u m 
is replaced by an integral: 

1 00 1 /*oo 

- ^ /(a, m )_+_ y du /( u ) for T^O, (27) 



m— — 00 
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and E then transforms into F by partial integration. 

It is convenient to put the position dependence of the dipoles in the atomic polarization function: a a (x) = 
~^2^a S(x — R) where R runs over all positions of the dipoles. For not including self-energy contributions, we 
have to replace D° by D' : 

, , . _ f D°(iw n , x) for x ^ . . 

M)(™>n,*)- j Q for f = Q • ^ 

The integration over x can be considered as matrix multiplications. Then the spectral densities of the energies (25,26) 
can be written as 

oo 1 

E(co) = -wd w ^2-MMr) (29) 

n=2 



FH = -jE^Tr((ao^D (30) 



The two energy densities are finite, since self-energy contributions (xt = x i+i ) are excluded by (28). 

Next we consider the limit to a continuous distribution of the dipoles. The sums over the positions are then replaced 
by integrals over the whole space. The polarization function is now 

a (iw n ,x) = a (iw n )9 K (x) , (31) 

where 0k{x) = 1 for a? in the region K occupied by the dielectric and otherwise. The polarization a (iw n ) is 
replaced by its density. D f is an integral operator with kernel D' (iui n ,x — y). The limit only exists if D f in (29, 
30) contains an ultraviolet cut-off in momentum space or a cut-off at small distances in x-space. This can be seen 
as follows. According to (29) (aoD' ) 2 must be a trace-class operator. This is true if a D' Q is a Hilbert - Schmidt 
operator, that means its integral kernel must be square integrable 

d 3 y J d 3 x\a (x)D' (x-y^\ 2 = \\a DXm=^-yol(K)\\DX,<oo. (32) 

Note that the square of a Hilbert - Schmidt operator is trace-class, and the same is true for all higher powers in (29, 
30). The dielectric must have a finite volume and the Green's function be square integrable. But a glance to (24) 
shows that an ultraviolet cut-off is necessary for this. The problem comes from the dipolc - dipole interaction. 

In x-space the dipole - dipole term $jj(aT) is not square integrable at x = 0. Therefore, we introduce a spatial 
cut-off at a minimum average distance a between the atoms as it was also done by Marachevsky [6]. We correct the 
Green's function D° for small distances as follows 

D' (lj,x) = D°(lu,x)9(x - a) = D°{uj, x) - G°{u, x), (33) 

where 

G°(w,x) = D°(oj,x)6(a-x). (34) 

We show next, that the splitting of D' = D° — G° leads to corresponding splitting in the energy E = En — Eq. In 
the sums (29,30) we introduce the first order term in a>o, which is since ao-Dd(O) = 0, and get 



- ]T -Tr((ao£>oD = Tr(ln(l - a ^° + aoG )) 

n—1 

= Tr(ln((l + a G°)(l - (1 + aoG )- 1 ^ ))) 

= Tr(ln(l+a G°)+ln(l-a£> )), (35) 
where a is the macroscopic polarizability given by 

a = (1 + a G r 1 ao- (36) 
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If we express ao by a in the first term of (35) we get 1 + a°G° = (1 — aG )^ 1 . We would like to split the trace in 
(35) into two separate traces. As discussed above this is only possible with an ultraviolet cut-off e~ Sk in the Green's 
functions. For the inner energy E we then get 

E = lim (4 - E 5 G ) (37) 
o — >u 

with 

1 °° 1 

E s D ( W )=wd w -J2-Tr((aD° 5 r) (38) 

n=l 

and 

E s G { W )=wdJ-f^±Tr{{aG°Y) (39) 

and similarly for the free energy F. In the limit S ^ E remains finite but E S D and Eq are both divergent. In the 
limit a — > also E is divergent. The energy E S D has the similar form as in (29), only that the atomic polarizability 
is replaced by the macroscopic polarizability a defined in (36) and D' is replaced by the cut-off Green's function D®. 
As the cut-off parameters a and 5 go to 0, we have 

G° 5 (u,x) ^ -±6(x) (40) 



and we get from (36) 



a(u,x,x') = ao } W >^ 5(x-x') (41) 
1 - ±a (u;,x) 



the Lorentz-Lorenz relation for the macroscopic polarizability. We shall use this relation in later calculations. 

E S D is for (5, a — > formally equivalent to the energy Ejj that is usually calculated in the Casimir problem. To bring 
the energy E S D in a convenient form, we perform the derivative in (38) using the relation 

wd w D° s = 2 (D° s + (D° s f) . (42) 

(38) then becomes 

E s D {w) = l -T^((wd w a) + 2(1 + a))D° (aD° s ) n , (43) 

n 

where the prime in the summation means that we just sum up in orders of a > 0. The last part ^2 n D°(aD°) n is the 
Green's function of the total electromagnetic field 

(Ei(g),Ej(x' )) =D\x-x') 

OO „ 

+ Y1 / d yi--- d ynDi kl (x-yi)a(y 1 )...a(y n )D° knj (y n -x'), (44) 

n=l J 

where a is given by (41). One can get this Green's function also by ordinary perturbation theory and Dyson resum- 
mation. It will be discussed extensively in the next section. Note that in (44) no cut-off is needed. 
With B = i/wV A E and the relation (sE(x) 2 ) = (B(x) 2 ) + div we finally get for energy density 

e D (w,x) = l(^(d w we(w,af))Ei(af)Ei(g)y + (B^B^f))') (45) 

where e = 1 + a and the prime means, that the 0-th order in a is suppressed. This corresponds only to the Too- 
component of electromagnetic stress tensor if there is no dispersion. The energy density in (45) was also used by 
Candelas [8] (cited by [9]). However this energy density is only well defined outside the dielectric. The total energy 
is infinite. 
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For two disconnected dielectrics it is possible to express the energy F by three parts: two containing just the energies 
of the isolated dielectrics and one containing the interaction energies of the to bodies. Let t»o(x) = <x${x) + a B (x), 
where A and B are the regions occupied by the two dielectrics. We then have 

F(w) = iTr(ln(l-K + ^K)) 

= l T r(ln(l - a$D' )) + ±Tr(ln(l - (1 - a£ D' )^ a B D> )) (46) 

The first term is just the energy Fa{w) = |Tr(ln(l — a A D' )^ of the dielectric in the region A. For the second term 
we use (1 - a A -Do) -1 = 1 + (1 - <* A D'q) -1 a A D' Q and we get 

F(w) = F A (w) + F B (w) + ^Tr(ln(l - (1 - a A D 1 ^ 1 a A D' Q {1 - a B D'^ 1 a B D',)) (47) 

where 

F B (w) = ±Tr(ln(l-a B D' )). (48) 

For the last expression in (47) we use 

(1 - a A (D° - G ))- 1 ^ = (1 - aAD^aA , (49) 

with ola = (1 — chqG^)^ 1 ^ . Since aAD'^as = ocaFPub for separated regions A and B, we get for the last term in 
(47) 

F AB {w) = ^(m(l - a A D A a B D B )) , (50) 

where 

D A = (l-D a A )- 1 D Q (51) 

and similar for Db- Da is just the Green's function of the electromagnetic field for only one dielectric in the region 
A as shall be shown in the next section. Hence we have for the free energy 

F(w) = F A (w) + F B (w) + Fab(w) (52) 

All the three energies are finite. For the energies Fa and F B one has to take the finite distance between the atoms 
into account. Fab is the only term where the distance between the dielectrics enters. Only this energy contributes to 
the Casimir force. It remains also finite in the limit a — > 0. Then for the polarizations a a and cub the Lorcntz-Lorenz 
formula (41) holds. 

To calculate the force between two dielectrics we will use the stress tensor of the electromagnetic field. The force is 
given by a contour integral over a region containing only one of the two bodies. Only the stress outside the dielectrics 
contribute to the force. It can be calculated as usual from the expression: 

T ik {x) - ^5 lk (E^x) E^x)}' - (E^x) E k (x)}' 

+ U lk (Bi (x) Bi(x))' - (Bi(x) B k (x))' (53) 

The prime again means that the 0-th order in a is suppressed. 



IV. BASIC INTEGRAL EQUATION 

To calculate the Casimir force from (53) one has to calculate the Green's function of the total electromagnetic field 
Dij(x, x') — (Ei(x) Ej(x')} given in (44). Let a(w, x) — a(w) 9k(x), where K is the region occupied by the dielectric. 
The Green's function can be expressed as an integral equation: 

Dij{u,x,x') = D°Jlj,x- x 1 ) + a(u>) / d 3 yD° k (u; 7 x - y)D kj (u,y,x"), (54) 

Jk 
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where no cut-off is necessary. This is our basic integral equation, ft separates the many-body problem from the 
geometry. The latter only enters in the region of integration K occupied by the dielectric. For simplicity we have 
assumed one kind of dielectrics and vacuum outside K, the generalization to more dielectrics is straightforward. 
From the integral equation one can derive a differential equation for D. From 

((cj 2 + A )6ik - did k )D kj (u, x) = -co 2 <%(x) (55) 

one easily derives 

((e(uj,x)lli 2 + A) Sik - did k ) D k j(u),x,x') = -uj 2 5 i:j {x - x' ) (56) 

The boundary conditions can also be derived from the integral equation (54). In D° only the part of the dipole-dipolc 
interaction is responsible for non-continuous boundary conditions. Hence the dipole-dipole approximation is necessary 
to produce hard boundary conditions and in the same time it is responsible for divergences for connected bodies, if 
the self-energy contributions arc ignored. 

A solution of the Casimir problem by solving the differential equation with boundary conditions was given by 
Lifshitz [7] for simple plain geometry. The solution of the boundary-value problem is much harder than the solution 
of (54). In principle the linear equation (54) can be solved numerically for any geometry. In the following we are 
interested in cases where this solution can be obtained analytically. We solve directly the integral equation without 
the use of any boundary condition. 

In the shorthand notation the basic integral equation writes as follows: 

D = D° + aD°9 K D (57) 
For a infinite dielectric (K = 1R 3 ) the equation can be solved in fc-space: 

D 1 = D° + aD°D 1 = {l-aD°)- 1 D° (58) 

leading to 



1 eLU 2 dn — kikj 1 

Let 



e R = \-e K (60) 

Substitute 9k in the basic integral equation (57) gives 

D = D a + aD°D-aD 6 K D. (61) 

Hence 

(1 - aD Q )D = D°- aD°9 K D (62) 

And using (58) we get 

D = D 1 — aD 1 6 K D (63) 
the complementary basic equation. To obtain the Casimir force we have to compute 

e K D0 K = 6 K D 1 6 K -a6 K D 1 e K D6 K . (64) 

Here we substitute (57) 

K D6 K = d K D l B K - a9 K D 1 6 K D e K - a 2 {9 K D 1 e K D )e K D6 K . (65) 

Now the same quantity as on the l.h.s has appeared, which allows the solution 

6 K D6 K = (1 + a 2 e K D 1 6 K D )- 1 9 K D 1 9 K (l - aD°9 K ). (66) 

The above formal solution (66) is very useful in the case where the operator in the inverse operates by simple 
multiplication: 

{9 K D 1 9 K D°)9 K D 1 9 K = 1 {w)9 K D 1 9 K . (67) 
As we will see this is indeed true for various simple geometries. Then we obtain from (66) and (57) the solutions 

9 K D9 K = / 9 K D 1 K (l-aD°0 K ), (68) 
1 + a z 7 

9 K D6 K = 9 R D°9 R + -^^9 R D°9 K D 1 9 R (l aD°9 R ). (69) 
1 + a A -y 
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V. PLAIN GEOMETRIES 

Let us first consider an infinitely thick plate in the region — oo < x < 0. In this case we choose the polarization 
vectors (11) as follows: Let n — (1,0,0), then 

ii(k) = p^- = -(0,-h,k 2 ) (70) 
\kAn\ P 

? /A fcA(fcAn) 1 . 2 , , , , . . 
|feA(feAn)| Pfc 

where p 2 = fc| + fcf and k 2 = k\+p 2 . D° (27) can be decomposed accordingly 

#° = £>? + £>§, (72) 

^^I dkl ^ lj W~4 eiklx ' (73) 

where S 2 , = w 2 — p 2 , w — iuj n , and 

D° 2ij (x) = 5{x)(SuZu - Sa) + ±jdki tebjr^-^e**. (74) 
The integrals can be easily computed: 

D° Ui {x) = && j % ^(e(x)e i >°* + 0(-x)e- is °*) (75) 

D° 2ij (x) = 5(x)(SuSij ~ Sij) + (6(x)e is °* + 6(-x) e - is °*) , (76) 

where £ 2 = ^/wp{p 2 ,id x k2,id x k^) due to the Fourier transformation. The corresponding formulae for the homoge- 
neous dielectric (59) are similar: 

Dhjix) = ^£^(5*)^ + 9{-x)e-™) (77) 

and 

j ? Fill 

D\ l0 (x) = -6{x)(Zutij - Sij) + (d(xy six + 6{-x)e-^ x ) , (78) 

where 

sl=ew 2 -p 2 , 6 ' = {p 2 , id x k 2 ,id x k 3 ). (79) 

The above decomposition into the two transverse modes (TE and TM mode) of the radiation field leads to a similar 
decomposition in (68) and (69). There is no mode mixing. Details of the calculation are given in Appendix A. For 
(68) we obtain 

6(-x)D(x,x')e(x') = + l ^ w * -f 2l &\6(-x)e-^e(x')e^\ (80) 

Vsi+so (si + so){p + siSq) ) 

and for the reflection (69) 

9(x)D(x,x')9(x') = 6{x)D°(x,x')6(x') 

+7^{ r i( w >P)b&j +r2(w, P )bi& j )0(x)e is ° x 9(x , )e is ° x ' (81) 
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where the second term herein represents the reflection at the boundary x — with the reflection factors r\ and r 2 for 
the TE and the TM mode given by 



and 



si-s ^e{w)w 2 -p 2 - yjw 2 -p 2 

ri(w,p) = = - (82) 

s i + s o y'e(w)w 2 — p 2 + \J w 2 — p 2 

s 1 -ss y/e(w)w 2 -p 2 - e{w)^/w 2 - p 2 

r 2 {w,p) = — — = — (83) 

si+eso y/e{w)w* — p z + e(w) x /w z — p l 

For later use we rename the Green's function for the region A : x < by D^a;, x'). Now we fill also the region 
B : x > a > with the dielectric with the same e, in between < ir < a we have vacuum. The Green's functions for 
the region B alone are obtained from the above formulae by substituting x by a — x: 

e B D°6 B = (^M, + (, Sl+ ^; 4 +5lS0 ) ^>(- - a)e-^0(a - *><-<•-'>, (84) 



9 s D(x, x')d s = 6{a - x)D°(x, x')9(a - x') 
iw 2 
2^ 



+ ^- (ri(«;,p)£iiffj + r 2 (w,p)Z 2l C 2j )0(a - x) e -^ x 8(a - x')e- is ° x ' e 2is ° a . 



(85) 



According to (57) the total Green's function is the solution of the integral equation 

D = D° + aD"9 A D + aD°6 B D. (86) 
Inserting the definitions of D A and D B , we get 

D = D A + aD A 6 B D (87) 
D = D B + aD B 6 A D. (88) 

Let us introduce the characteristic function 

6 C = 1 - 9 A - 6 B (89) 
for the region between the two dielectrics. Then by combining the equations we find for the polarization A: 

e A D x e c = e A D A e c + a.e A D A e B D B e c + a 2 e A D A e B D B e A D x e c 

= - — 1 -^{6 A D A e c + ae A D A e B D B e c ). (90) 

Here we have used the fact that the operator 

(6 A D A 6 B D B )6 A D A 6 B = jx6aD a 6 b (91) 
operates simply by multiplication, this is shown in Appendix A. Substitution into (88) gives 

e c D x 8 c = e c D B c + n x a (6 c D B 6 A D A 6c + a6 c D B 9 A D A 6 B D B e c ) , (92) 

where 

1 

1 - a 2 -f X ' 



MA = q — • (93) 



Substitution of D B by its defining equation D B = D x + aD x 9 B D B leads to 

An , .. „a ndfl nB( 



C D x e c = 6 c D x e c + n x ad c D x e A D A 9 c + Vxa9 c D° x e B D B d c 

+ i2xa 2 e c D° x 9 A D A 9 B D B 9 c + ^ xa 2 e c D x e B D B 6 A D A 6 c (94) 
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The calculation of (94) is given in Appendix A with the following results: 

D x (x,x') = D x (x,x') 

+ Ma r x ^ (e is ^+^ + e -»o(x+x'-2a)^ 

2 

+ ma r\ e 2lSoa ^- £xCx (e iso(x - x,) + e -<»o(*-*') , j , (95) 



^ = !_ r 2 e 2 tSo a ^ ( 96 ) 

with A = 1, 2 for the two modes and where r\ are the reflection factors given in (82) and (83). Since the 0-th order 
does not contribute later on, we introduce D w = D — D°: 

D%(w; x, x') = duix^ix'yfiw; x, x') + ^i{x)^j {x')*f 2 {w; x, x'), (97) 

where 

£i(k) = -(0,-k 3 ,k 2 ) (98) 

= — (p 2 ,id x k 2 ,id x k 3 ). (99) 
wp 

VI. CASIMIR FORCE FOR PLAIN GEOMETRY 

Let us assume for a moment that the dielectric A covers the half-plain — oo < x < and dielectric B the region 
a < x < R. Then the force per unit square on the surface of B is equal to the difference of the stress tensors 
Tu(a—) — Tu(R+). The term Tn(R+) vanishes, but Tu(a— ) contains contributions from the reflection at x = R. 
These contributions vanish in the limit R — > oo. Then only the Green's function for arguments between the two 
dielectrics contributes, which we have calculated in the last section. This argument is necessary because otherwise 
the stress tensor at infinity gives rise to an infinite "volume force" . 
The stress tensor for 1-direction is given by 

Tn(f) = \{E 2 {xf + E 3 (xf - E^xf 

+ B 2 {xf + B 3 (x) 2 - B 1 (x) 2 )' . (100) 
To find the Green's functions for the magnetic field, we use the field equations. From 

E = -d t A and B = V A A (101) 



we get for the Fourier transforms 



By (70) (71) we have 



B(k) = —k A E(k). (102) 



-k A = &(k) and -kA&(k) = -£i(k), (103) 



using k 2 = uj 2 . Now we are ready to calculate the Casimir force. We obtain with (97) 

Tll(f) = ^£ e " mT (27)2 / ^^ 2 + 9A / - P 2 )(f 1 (w;x,x') + f 2 (w;x,x')) 



(104) 
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where the fj are defined in (97). Only the terms depending on the difference (x — x') in fj contribute because they 
are proportional to 

w 2 + si -p 2 = 2sl (105) 

This leads to the following final result 

where n and ri are the reflection factors for the two modes given by (82) and (83). The result (106) is the Lifshitz 
formula [4] [7]. 

VII. SPHERICAL GEOMETRY 

As mentioned in Sect. 3 for compact dielectrics the energy Ed is always divergent, since it ignores the self-energy 
contribution Eg in equation (37). Nevertheless the above method to calculate the photon Green's function still works 
for spherical geometries as shall be shown in this section. As a consequence, it is not hard to compute the Casimir 
force on the spherical shell and rederive the result of ref. [10]. 

The region occupied by the dielectric is now the sphere 

K = {x etf |x 2 < R 2 }. (107) 
The free Green's function (23) (24) is then written in spherical coordinates (r, <p) using the expansion 

r k 

Lm 

where ji are spherical Bessel functions and YJ m the spherical harmonics. For f /x' we then have 

OG 

D° f (w,x,x>) = - £ ^(f)^(f') / dk-z^^MkrMWYr&tfYr&rfy ■ (109) 

\=1,2 q lm 

The fc-integral can be extended to (— oo, +oo) and then computed using the theorem of residues, taking w — iw n (22) 
into account. The result is 



iw 3 [6(r - r')hf ] \wr)jt(wr') + 6{r' - r)ji(wr)hf ] \wr'j\ , (110) 



where ftj are the Bessel functions of the third kind. The two transversal polarization vectors perpendicular to 
k <~ p= —id/dx are constructed as follows: First we choose 

fi(x) = — 1 xhp= , l :L. (Ill) 
This is the angular momentum operator which gives the so-called electric TE mode; it is normalized according to 



The second polarization vector is then 



6(x)= JL= pA£. (113) 

Wy/l[l + 1) 
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This is the TM mode which is also normalized because |fc| = w. Since the metric tensor in spherical coordinates 
is non-trivial, it is convenient to distinguish between upper and lower tensor components. Let Q — 1(1 + 1). In 
components we have 

a(r,4>,e) = -L_L_(o,-a 9) ^) 

y/Q r sin ft V / 
Su(r,<f>,0) = -^r(o, -sin#d e , (114) 

ky/Q V r r 2 sin r z J 
£,2i{r,(j>,0) = -j-^ , d r r 8$ , d r r d e ^J (115) 

We have orthogonality in the sense [10] 

00 

/ (&(mr)Yr(*))\ii(Z)9(r)Y l r\*)d 3 x = S w 5 mm> J drr 2 .r(r)g(r) (116) 



&Wf{r)Yr{-)) ^ l (S) 9 (r)Y^\-)d 3 x^du'6 mm> — ldr[-(rn^(r9)+l(l + l)rg , (117) 



where we sum over the lower and upper indices i. The mixed expressions are 0. This has the important consequence 
that different Vs and polarizations don't mix. Therefore, only the radial integrations are non-trivial. The free Green's 
function (109) is now equal to 



D°/(w,x,x") = ]T ^mi{x')Y,Yr^MYi m ^^rD {r,r'), 



(118) 



A=l,2 



Irn 



(119) 



(120) 



where 

D°(r,r') = iw a 9(r - r')hf ] ' (wr)jt(wr') + 6(r' - r)/i ; (1) (wr')ji(wr) 
and for the infinite dielectric we have 

D 1 (w, x, x') = -D°(y/ew,x,x'). 

The essential point in the solution of the problem is the calculation of (67). For the TE mode and r, r' < R we get 

00 

e K D 1 e k D°e K \ l {r,r') = J dg g 2 D 1 (r, g)D°(g,r') 

R 

= -(h\ 1 \h\ 1) )i^°ji(Viwr)j l {wr f ) , 



(121) 



where we have introduced the non-symmetric bilinear form of Bessel functions: 

R 



(fl,9l)i = J dgg 2 fi(wg)gi(y/ewg) 



R 



aw' 



wf' l (wR)gi(<JewR) - y /ewfi(wR)g' l (s/ewR) 



Then (67) operates multiplicatively, indeed, 

{e K D 1 e k D )e K D 1 ef i \ l = iu0kD 1 6k\ v 



(122) 



(123) 
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where 

7i/ = ^ 6 W 1) .M 1) )i0^i«)i- ( 12 4) 

Since 

(9 R D 1 9 R -a9 K D 1 9 R D°9 R )\ l = -a^w 6 ■(j l ,h\ 1) ) 1 j l (y/i W r)h < i 1) \wr') , (125) 
(see also (All) in the appendix) we obtain for (68) 

0KDfe R \ x=1 = -—^-ay/hv^h^h 



for A = 1 and similarly for A = 2 

V**L=2 - l + a 2 72 "^ 



x ^.•(^'(fo^w^r^'^OiK^)^^'), (127) 



with 



-2 = ^0'l,Ji>2<M 1) .M 1) >2. 



where we have introduced a second bilinear form for Bessel functions 

R 



(f i, 91)2 = I dr ^(rfi(wr)^^(rgi(^wr)^+l(l + l)fi(wr)gi(y/ewr) 







d 



d ^(RMwR))]Rgi(V^wR) + w 2 (f l ,g l ) 1 . (129) 
For (69) we obtain: 



lm 



where for A = 1 



a 2 Vsw 6 (j l Ji} 1 (ji,h ( l 1) } 1 
fJ-il = Tn m — ' ( ic5i J 

and similarly for A = 2: 

m = — 2 TT-, 7TT— • 132) 

l + O 2 ^0j,J«>2<fc I (1) ,M ) >2 

VIII. COMMENTS ON THE METALLIC LIMIT 

In this section the metallic limit e — > is considered. By the Lorentz-Lorenz relation (41) this means to put the 
atomic polarization a — 3. It is not clear whether in our microscopic theory the perturbation expansion for the 
Green's function converges for ao = 3. Our microscopic theory does not give an adequate description for metals. In 
fact we'll show below that in infinite homogeneous dielectrics the series only converges for a < §. Nevertheless for 
the Casimir effect only Green's function that crosses the boundaries contributes to the force. Indeed the series then 
converges in the TE mode, at least, for small frequencies p 2 > 2u 2 , where u = w n = —iw. We'll show this only for 
the special case of fiat geometries. 
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The Green's function of the radiation field is given by 

D = D° + a D a 9 A D° + a 2 D°9 A D'9 A D , (133) 

where 

D' = D' + a D' 6 A D' + a 2 (D' 6 A ) 2 £>£ + ... (134) 
with D' a = D° + i. In momentum-space D' is given by 

r\l r ■ r\ ~~U 2 6ij — kikj 1 

1 2ii "I - Z' c k^kj \ 2 k%kj 



3 fc 2 +« 2 * 2 J 3 fc 2 (135) 

It is clear that in infinite homogeneous dielectrics the perturbation expansion only converges for < | because of 
the dipolc-dipolc interaction. 

For the Green's function that crosses the boundaries, let us first consider the case u = 0. Then only the dipole-dipolc 
interaction contributes. We have 

D o(0, fc) = + ^% = D l (k) + Isy (136) 

Let A be the region of the right half-plain with < x < +oo. The longitudinal Green's function is given by 

D l (x -y) = -(-id x ,p 2 ,p 3 ) T (id y ,p 2 ,p 3 )d°(x - y) , (137) 

where 

d ° {x) = hj ^IzjTpt^ 1 * = Yp (°^ PX + 0(-x)e px ) , (138) 



with p — \Jp\ + p 2 - Then a straight-forward calculation gives 



9 A D l 9 A D l 9 A = -±e A D l Q 6 A (139) 



For the Green's function 9 A D9 A we get 



This is convergent for ao < 6. Hence we get 



^ = (i-f(i-f + Qf + ...))e A Dle A (mo) 



(9 A I>^ = (1 - ^ T -^.)(9 A £>^ (141) 

This is for a — 3 as it should for metals. For the reflection we get with (^6 A D 6 A D 6 A )(x, x') = — ^D l (x + x') for 
x, x' > 

D(x, x') =D l (x- x') ^-^L-D^x + x') =D l (x- x') ^D^x + x') (142) 

where 

e = 1 + T^k- (143) 

1 3 

For u > it is much harder to show that series converges in the metallic limit. We will show this only for the 
special case where p 2 > 2u 2 > and only for TE mode. Let the atomic polarization be any function with 

a (iu)<3 for m > , a (0) = 3. (144) 
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We express the expansion for the Green's function D' (134) formally by the integral equation 

D' = D' + a D' 6 A D' , (145) 
where the TE mode of the free Green's function D' is given by 

D' {iu,k)=D\u,k) + \ = v ^+ 1 -. (146) 

For p 2 = k\ + k\ > 2u 2 > the expansion for Green's function D[ in a infinite homogeneous dielectric converges and 
is given by 

a a 2 



with 



D[(iu,k) = — + "£)!(«, fc), (147) 



-ii 2 

D\iu,k) = WT ^ 2 . (148) 



This leads to the complementary integral equation 

D' = D[ - a D[e A D' . (149) 

We then get 



6 A D'e A = 8 A D 9 A + a 8 A D°9 A D'9 A + ^-9 A D'9 A 



= 9 A D°9 A + a Q 9 A D"D'9 A - a 9 A D°9 A D'9 A + -^9 A D'9 A . (150) 

Substitution of D' (149) into the second term of the second line and use of D° + a a D°D' 1 = -^D 1 gives 

9 A D'9 A = — 9 A D 1 8 A -a8 A D 1 9 A D l 9 A + ^-9 A D'9 A 
a 3 

= — 9 A D X B A - a9 A D 1 9 A D°9 A - aoa9 A D 1 6 A D°6 A D'9 A + ^9 A D'9 A (151) 
a 3 

Now the operator 

(-a a9 A D 1 9 A D a 9 A + ^-)9 A D 1 6 A - 1 8 A D X 9 A (152) 

operates by a simple multiplication. From (A9) in the appendix we get 

go_°o('i-*o) 2 <! (153) 
3 a4s\so 

for u > 0, hence the perturbation expansion converges. The reflection factor for the TE mode is then given by (82) 



/ \ s-l-sq ^eu 2 +p 2 - ^Ju 2 +p 2 

r 1 (u,p) = — = —*_=_ ■ (154) 

si + s yjsu 2 + p 2 + V?i 2 + p 2 

This formula for the reflection for the TE mode only holds for p 2 > 2u 2 > 0. For u = the TE mode does not 
contribute to the reflection and to the free energy independent of the model adopted for the dielectric function. 

This cannot be deduced directly from the Lifshitz formula, where the TE mode does not contribute only if 
lim u ^ e{iu)u 2 — 0, which does not hold in the plasma model for metals where 

u 2 u 2 

ao(iu) = ^ — 2 ► e{iu) = H — | . (155) 

u 2 + ^ u 

Although lim u _^ r\{u,p) ^ for the plasma model, in our microscopic theory the TE mode does not contribute to 
the free energy for u = 0. 

A different result is also obtained for the TM mode for the reflection at a finite metallic plate in the region 
> x > —R. It can be seen from formula (155), that for u — > the atomic polarization ao(0) = 3 does not depend 
from the plasma frequency u p . This has the consequence that as can be shown, that the reflection at the plate does 
not depend on the thickness R of the plate for u = and we get formula (142) with alpha = 3. Whereas in the 
macroscopic theory there is a contribution from the reflection at x = R. 

The small frequency behavior is important for the temperature dependence of the Casimir force, that is widely 
discussed in recent papers [11-17]. 
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APPENDIX A: 



In this Appendix wc do the main calculations for the plain geometries. All the calculations are simple since it 
suffices to calculate surface integrals. First we have to compute the reflection at one plate. We use formula (66) 

6 K D6 R = (1 + a^KD^nD^-^KD^jiil - aD a 9 R ). (Al) 

To compute the integral 

' d 3 y9 K (x)D 1 (x-y)D°(y-z)9 K (z) (A2) 



we use the differential equations 



(A +u 2 )D%{y) = -(didj+rfSMff) 
(A +ew 2 )D\ j {y) = -1(9^- W<^W)- (A3) 



On the r.h.s is a function with point support. Since in (A2) y and x resp. z have different support D° and D 1 simply 
solve the wave equations. By Green's theorem we get a surface integral 

/ d 3 ye K {S)D\x-y)D\y- z)9 K (z) = -L f d ^(y) 9 K (x) 

Jit auJ JdK 

x ( y (V^D\x-y))D%y-z^-D\x-y)V g D\y-z^e K (z) (A4) 

where the surface vector <j(y) shows inside the region K. Because of the hight symmetry of the plain or the spherical 

geometries no integral at all has to be computed. 

For plain geometries and A representing the region x < we get 

9 A B\9 A B\9 A = lUj2(s l - g^ &fi M-x)e- is ^6{-x')e- is ° x ' (A5) 
for the electric TE mode and 

e A D\e A Dle A = luj2{si ~ So)(p ' zJi^ ^.e { - x )e-^e{-x')e-^' (A6) 



for the magnetic TM mode. Similarly we get 



e A D\e A Dle A D\e A = {s \ So) x ^ u ge(-x)6(x')e- is ^-^ (at) 

for the TE mode and 

e A D\e A Dle A Dle A = (fi^lVzW x ^ 2i e 2j e(-x)e(x')e-^-^ (as) 

creu;* 4siSo 2si J 

for the TM mode. Hence 9 A D 1 9 A D°9 A operates by a simple multiplication on 9 A D 1 9 A with the factors 

Oi - s o) 2 , (si - s ) 2 (p 2 - sis ) 2 
7i = Hi and 72 = — — (A9 

for the two modes. Now we compute formula (68) 

9 A D9 A = —L^-(9 A D 1 9 A - a>9 A D 1 9 A D°9 A ). (A10) 

The expression in the bracket can again be computed by the use of the differential equations (A3). The Dirac functions 
in the differential equations eliminates the first part. With Green's theorem we obtain 

9 A D l 9 A - a9 A D 1 9 A D"9 A 

= / d 3 y6 A (x) (A yD 1 (x — y) D°(y — x') — D x {x — y) A y D°(y — x')) 6 A (x') 

J A 

= -\t ^s{g)e A {^(^ y D\s-^{i}-s')-D\3-f^ y D ^-s'))e A {s') (ah) 
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For plain geometry we obtain 

e A D\e A {\ - r%o A ) = tuj2 [ Sl + s ^ u eiA- x ) e ~ islX() ( x ') eiS0X ' : A 1 -> 

and 



e A D\e A (i Dle A ) = tuj2{s \ + So)i f ± SlS2 \ 2l c 2 A^ lSlX0 ^'y saX ' 



(A13) 



The factors in (A10) simplifies to 



f 4 Sl s 1 4s lSo ew 4 

aIld 1 , „,2- = /„ , I X2^2 I „ „ X2 ■ ( AU ) 



1 + a 2 7i (si + s ) 2 ' 1 + a 2 72 (si + s ) 2 (p 2 + sis ) 2 
where in the second factor ecu 4 = (s\ + p 2 )(s^ + p 2 ) was used in the calculation. Hence we get for (Al) 

2 ■ / — 4 

e(-x)D A (x,x')6(x') = f-^—^l + W " W -i 2l Q 3 )9{-x)e-^e{x')e^' (A15) 

Vsi + so (si + so)(p + siso) > 

For the reflection at one plate we have to compute (69): 

9 A D9 A = 9 A D°9 A + a 9 A D°9 A D 1 9 A (l - aD°9 A ) . (A16) 

1 + CT7 

The further integration adds the factors -^gf 1 and - (si ^° o ) ^!j 2 SlS ° ) . We obtain the result (81) for the reflection at 
one plate: 

2 

6(x)D A (x,x')e(x') = 9(x)D°(x,x')e(x') + ^- (n(tu, + r2{w,p)^>j6{x)e is ° x e{x l )e is ^' (A17) 

with the two reflection factors 

Si-sq (si - s )(p 2 - SiS ) , . 

n = and r 2 = -7 r . (A18) 

si+s (si + s ){p 2 + S1S0) 

The reflection factor r 2 for the TM mode simplifies with s\ = ew 2 — p 2 and s 2 , = w 2 — p 2 to 

r 2 = ■ (A19) 

si + es 

For a plate in the region B with a; > a we just have to replace in the formulae x by a — x. We will just need the 
formula for 9 B D B 9g. From (A15) we get 

6(x-a)D B (x,x')e(a-x') = (-^-^ + rfr&W - a)e^- Q )% - a/)e is ° . (A20) 

For the Green's function for two parallel plates we use formula (94) 

9 C D X 9 C = 0cD o x ec + »xa6 c Dle A D A 6c + Vxa6cD x 9 B D B 9c 

+ i2 X a 2 9 c D° x 9 A D A 9 B D B 9 c + i2 X a 2 9 c D° x 9 B D B 9 A D A 9 c (A21) 

where 

^ = 1 - 2 ■ (A22) 
1 - a 2 -f X 

and 7a is given by 

(9 A D A 9 B D B )9 A D A 9 B = lx 9 A D A 9 B (A23) 
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First we compute the expression in the quote in (A23). For D A and D B similar differential equation holds as in (A3), 
so that we just have to compute a surface integral giving 

e A D^e B Dfe A = e «»o° ^* 1 e(-x) e - lsix e(-x')e~ lsox ' (A24) 

for the TE mode and 

e A D£e B D*e A = zb/p^i L: liffi > e Q '°° 6^(-z)e- islx 0(-z')e- <so3: ' (A25) 

a(si + so) (P + siso) 

for the TM mode. A further integration adds the factors ' ( ^J 2 So) and ^"^^r^ , so that we get for 7A in (A23) 

a 2 7! = ^-^e— and = fc^V^jifgg ^. (A2g) 

(si+so) 2 (si + sor(p 2 + siso)- 4 

where the reflection factors from (A18) appear again. For /Ua we get 

^ = 1 _ (A27) 

Finally we have to compute (A21). The second term is just the reflection at the plate A and is given (A17). The 
third term is the reflection at the plate B, where x has to be substituted by a — x. For the forth and fifth terms the 

further integration adds again the factors and t(si ~ s °%: SlSa) to (A24) and (A25). This gives for x and x' 

between the two plates 



D x (x,x') = D a x (x-x') 

2 

+ A* A r \ ( e »o(*+*') + e i S0 (2a-x-x')^ 

■ 2 

+ A* a r 2 x e 2lSoa l ^-txC\ (e lSa{x ~ x>) + e^'^'A (A28) 

2sn V / 



for A = 1, 2 and with ^ A given in (A27) and the reflection factors r\ given in (A18). 
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